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ABSTRACT 

The  results  of  three  quasi-objective  techniques  for 
enhancing  the  accuracy  of  24-hour  cloud  cover  forecasts 
in  sparse  data  areas, using  satellite  data  and  numerical- 
ly-produced analyses  and  prognoses,  are  presented.   All 
three  of  the  techniques  involve  the  modification  of  the 
analyzed  500-mb  relative  vorticity,  as  based  on  the  sat- 
ellite-observed cloud  patterns   and  the  associated 
changes  to  the  prognostic  fields  of  500-mb  relative  vor- 
ticity.  The  final  24-hour  cloud-cover  forecast  is  based 
on  the  manually-modified  relative  vorticity  pattern   and 
three  different  advecting  currents.   These  advecting 
currents  are  based  on  the  prognostic  500-mb  height  field, 
modified  and  unmodified,  and  the  prognostic  500-mb  space 
mean  field.   The  study  was  conducted  over  a  large  area 
of  the  North  Pacific  Ocean  for  the  period  9  March  to 
13  April  1972. 
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I.   INTRODUCTION 


The  objective  of  this  study  was  to  investigate  the 
possibility  of  increasing  the  accuracy  of  existing  cloud 
cover  forecasting  techniques  in  sparse-data  areas  by  us- 
ing satellite  data  in  conjunction  with  numerical  analyses 
and  prognoses.   In  particular,  an  improvement  was  sought 
on  Bittner's  [1]  quasi-objective  method  of  using  the  ad- 
vection  of  relative  vorticity  on  the  500-mb  surface  to 
determine  the  primary  areas  of  cloud  cover. 

Bittner's  quasi-objective  method  uses  current  satel- 
lite-observed cloud  patterns  to  modify  the  analyzed  pat- 
terns of  relative  vorticity.   The  difference  between  the 
initial  and  modified  relative  vorticity  patterns  is  then 
subjectively  applied  to  the  24-  or  36-hour  prognostic 
patterns  of  relative  vorticity.   The  primary  areas  of 
cloud  cover  are  then  related  to  the  areas  of  positive 
relative  vorticity  advection  (PVA)  as  determined  by  the 
unmodified  24-  or  36-hour  prosnostic  500-mb  height 
fields  (Method  1) . 

Theoretically,  an  improvement  in  the  cloud  cover  fore- 
casts could  be  made  if  the  prognostic  500-mb  height 
fields  were  modified  by  the  amounts  associated  with  the 
changes  in  the  prognostic  relative  vorticity  fields  prior 
to  determining  the  areas  of  positive  relative  vorticity 
advection  (Method  2). 


Also,  the  theoretical  development  of  Fjortoft  [2]  In 
manually  producing  a  barotropic  24-hour  forecast  of  the 
500-mb  height  values  indicates  that  using  prognostic 
space-mean  fields  to  determine  the  areas  of  positive 
relative  vorticity  advection  might  show  an  improvement 
in  the  verification  of  24-hour  cloud  cover  forecasts 
(Method  3)  . 

Therefore,  this  study  was  specifically  undertaken  to 
determine  which  of  the  above  three  methods  are  the  most 
accurate  for  forecasting  24-hour  synoptic-scaled  cloud 
patterns.   All  three  techniques  were  used  simultaneously 
to  produce  cloud-cover  forecasts  over  the  North  Pacific 
Ocean  from  20°N  to  50°N  and  from  170°E  to  125°W.   The 
forecast  period  extended  from  9  March  to  13  April  1972. 

Modifying  the  500-mb  height  prognoses  in  addition  to 
the  relative  vorticity  prognoses  took  slightly  longer 
than  the  other  two  methods,  but  the  resultant  improvement 
in  the  forecasts  seemed  to  indicate  that  the  small  amount 
of  additional  work  was  well  worth  the  effort  involved. 


II.   BACKGROUND 

The  oceanic  areas  of  the  world  are  relatively  void 
of  conventional  meteorological  observations,  especially 
above  the  earth's  surface.   Over  the  middle-latitude 
areas  the  synoptic-scaled  cloud  patterns  attain  their 
best  development.   Fortunately,  these  cloud  patterns, 
which  are  so  revealing  of  synoptic  systems  to  trained 
meteorologists,  are  easily  observed  with  satellites. 
Consequently,  for  the  past  several  years,  meteorologists 
have  striven  to  incorporate  the  vast  quantity  of  synoptic- 
scaled  information  obtained  from  satellites  into  useful 
models  for  improved  forecasts.   Further,  there  has  been 
considerable  effort  to  incorporate  this  satellite-derived 
data  into  numerical  analysis  schemes,  as  Nagle  and  Hayden[7], 
Nagle  and  Clark  [6],  Mantei  and  Workman  [5],  Zeigler  (8] 
and  Glaes  [3].   Recently  Bittner  [1]  has  illustrated  a 
method  of  using  these  satellite-observed  synoptic-scaled 
cloud  patterns  to  manually  improve  the  numerical  analysis 
of  the  500-mb  relative  vorticity  field.   He  then  carried 
this  improvement  forward  24  to  36  hours  by  manually  adjust- 
ing the  numerically-produced  prognostic  500-mb  relative 
vorticity  field.   Finally,  he  converted  this  modified  prog- 
nostic field  into  cloud  and  precipitation  forecasts. 

In  an  equivalent  barotropic  numerical  model  of  the 
atmosphere  it  may  be  easily  shown  that  the  sign  of  the 
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advection  of  relative  vorticity  at  the  level  of  nondiver- 
gence  (LND)  is  the  same  as  the  sign  of  the  vertical  motion 
through  this  level.   In  fact,  in  this  simple  model  of  the 
atmosphere  the  sign  of  the  vertical  motion  is  the  same 
throughout  a  given  column  of  air,  if  terrain  effects  are 
neglected.   Therefore,  positive  advection  of  relative  vor- 
ticity at  the  LND  produces  positive  (upward)  vertical 
motion  throughout  the  column  of  air.   In  the  real  atmosphere 
the  sign  of  the  advection  of  relative  vorticity,  or  geo- 
strophic  relative  vorticity,  on  the  500-mb  pressure  surface 
is  highly  correlated  with  the  sign  of  the  synoptic-scaled 
vertical  motion  through  this  surface.   That  is,  if  the 
advection  of  geostrophic  relative  vorticity  on  the  500-mb 
surface  is  negative  then  the  synoptic-scaled  vertical  motion 
in  the  vicinity  of  the  500-mb  level  is  also  negative  (i.e. 
subsidence)  almost  all  of  the  time. 

If  hygroscopic  nuclei  and  moisture  are  present  within 
the  atmosphere  then  the  formation  and  dissipation  of 
clouds  are  largely  determined  by  the  sign  of  the  vertical 
motion.   That  is,  with  the  above  conditions,  positive 
vertical  motion  on  a  synoptic  scale  will  produce  synoptic- 
scaled  cloud  patterns  when  the  time  scale  is  sufficient 
to  cover  the  phase  lag  between  initial  lifting  and  subse- 
quent condensation.   With  a  given  moisture  distribution, 
the  stronger  the  advection  of  vorticity  the  shorter  the 
time  lag  between  the  vertical  ascent  and  its  tell-tale 
signature,  the  synoptic-scaled  cloud  patterns. 
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For  a  given  analysis  or  prognosis  of  the  height  field 
of  a  constant  pressure  surface  the  relative  vorticity  may 
be  obtained  diagnos t ic ally  from  an  equation  which  relates 
the  mass  field  to  the  wind  field.   The  exact  form  of  this 
diagnostic  equation  depends  upon  the  assumptions  made 
concerning  the  relationship  between  the  wind  field  and 
the  mass  field.   Fjortoft  [2]  demonstrated  that  the  height 
field  of  a  constant  pressure  surface  may  be  smoothed  by 
using  a  space  differential  operator  called  the  Laplacian 
or  del-squared  operator.   This  smoothed  height  field  is 
more  commonly  known  in  meteorology  as  the  space  mean.   The 
space  mean  represents  the  initial  height  field  after  the 
shorter  wave  lengths  have  been  selectively  removed.   The 
exact  short  wave  length  which  is  completely  removed  from 
this  field  depends  upon  the  grid  length  selected  for  the 
finite  difference  equation  used  in  the  approximation  of 
the  analytical  equation.   Fjortoft  proved  that,  except  for 
a  latitude-dependent  weighting  function,  the  geostrophic 
relative  vorticity  field  is  obtained  when  the  initial 
height  field  is  subtracted  from  this  space-mean  height 
field. 

The  height  field  of  a  constant  pressure  surface  which 
is  obtained  through  numerical  analysis  or  prognosis  may 
be  subjected  to  the  scale  separation  technique  of  Holl  [4] 
In  effect  the  initial  height  field  is  repeatedly  smoothed 
with  a  Laplacian  operator  until  a  residual  field  called 
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the  SR  field  remains.   This  SR  field  may  be  roughly  con- 
strued to  be  equivalent  to  the  space-mean  field  described 
above.   When  this  SR  field  is  subtracted  from  the  initial 
height  only  a  disturbance  field,  called  an  SD  field,  re- 
mains.  The  SD  field  roughly  represents  the  combined  effect 
of  all  the  numerically  repr esentable  short  waves  with  wave 
lengths  shorter  than  wave  number  six.   Likewise  the  SR  pat- 
tern approximately  depicts  wave  numbers  zero  (polar  vortex) 
through  six.   Note  that  the  field  which  approximates  the 
geostrophic  relative  vorticity  is  the  space-mean  minus  the 
initial  height  field  while  the  SD  field  is  minus  the  SR 
(space  mean)  plus  the  initial  height  field.   That  is,  when 
the  pattern  of  the  SD  analysis  is  used  to  approximate  the 
pattern  of  the  geostrophic  relative  vorticity,  the  sign  of 
the  SD  centers  must  be  reversed.   See  Appendix  A. 

Consequently,  the  above  concepts  led  to  the  development 
of  quasi-objective  rules  for  relating  synoptic-scale  clouds 
to  the  500-mb  SD  field.   Bittner  [1]  outlined  these  and  for 
convenience  they  are  listed  in  Appendix  B.   Bittner's 
method  consists  of  using  the  most  recent  satellite-observed 
cloud  patterns  to  modify  the  current  500-mb  SD  analysis. 
The  resulting  difference  field  between  the  initial  and  modi- 
fied SD  patterns  is  applied  to  the  prognostic  SD  field  to 
produce  modified  forecast  SD  patterns. 

The  particular  step  of  modifying  the  SD  patterns  of  the 
prognostic  field  based  on  the  modifications  applied  to  the 
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current  SD  patterns  appears  to  be  the  most  subjective  of 
all  the  procedural  steps  of  Bittner's  method.   If  there 
is  no  energy  transformation  occurring  between  the  short 
waves  represented  by  the  SD  field  and  the  long  waves  por- 
trayed by  the  SR  field  during  the  period  of  the  prognosis, 
then  the  modifications  of  the  prognostic  SD  patterns  are 
rather  straight  forward.   However,  for  a  given  cyclonic 
or  anticyclonic  system  within  an  initial  SD  field,  the 
greater  the  forecast  change  of  development  of  a  short-wave 
feature  the  less  reliable  the  modification  which  is  applied 
to  the  SD  system  of  the  prognostic  field.   That  is,  the 
more  an  SD  center  changes  during  a  prognostic  period,  the 
more  difficult  it  becomes  to  propagate  an  initial  modifi- 
cation of  the  SD  center  to  its  prognostic  counterpart. 

The  modified  SD  prognosis  is  then  used  in  conjunction 
with  the  unchanged  500-mb  height  prognosis  to  derive  a 
subjective  cloud  and  precipitation  forecast  valid  at  the 
time  of  the  numerical  height  forecast.   Bittner  makes  the 
following  assumptions  relative  to  his  synoptic-scaled 
cloud  pattern  forecasts: 

1)  Large-scale  vertical  motion  fields  are  highly  cor- 
related with  the  characteristic  distribution  of  extratropi- 
cal,  synoptic-scale  cloud  patterns  as  seen  by  satellites. 

2)  In  turn,  the  vertical  motion  patterns  are  closely 
related  to  the  advection  of  relative  vorticity  at  500  mb . 

3)  An  SD  field  is  equivalent  to  a  relative  vorticity 
field.   See  Appendix  A. 
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A)   Any  errors  in  the  SD  field  based  on  the  initial 
500-mb  analysis  are  propagated  during  the  numerical  prog- 
nosis to  produce  similar  errors  in  the  SD  field  based  on 
the  prognostic  500-mb  heights.   Hence  the  error  field  is 
assumed  to  be  conserved  during  the  prognosis. 
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III.   DATA  SOURCES  AND  SELECTION  OF  STUDY  AREA  AND  PERIOD 

An  area  of  the  North  Pacific  Ocean  extending  from  20°N 
to  50°N  and  from  170°E  to  125°W  was  selected  for  investiga- 
tion for  several  reasons. 

This  area  was  selected  primarily  because  it  is  both 
analyzed  numerically  by  FNWC  and  has  almost  continuous  day- 
light picture  coverage  by  the  ATS-1  satellite.   ATS-1  is 
the  designator  for  the  geostationary  Applications  Technology 
Satellite  positioned  near  the  equator  of  the  Pacific  Ocean 
since  November  1967.   Its  spin-scan  cloud  camera  takes  time- 
lapse  photographs  approximately  every  25  minutes.   The  re- 
sulting pictures  cover  a  region  from  about  50°S  to  50°N  over 
a  span  of  100  degrees  of  longitude  centered  at  150°W,  with 
a  resolution  of  two  miles  at  the  center. 

The  subject  ATS-1  photographs  were  generally  available 
locally  from  FNWC  on  a  real-time  basis  Monday  through  Friday 
Since  the  Saturday  photograph  was  needed  for  verification 
of  Friday's  forecast,  the  study  period  was  usually  limited 
to  four  days  per  week,  namely  Monday  through  Thursday.  Some 
of  these  photographs  and  others  that  were  missing  on  occa- 
sion were  supplied  by  the  Naval  Air  Systems  Command's  Pro- 
ject FAMOS.    The  FNWC  computer  products  were  also  available 
locally  on  a  real-time  basis. 


Presently  part  of  the  Environmental  Prediction  Research 
Facility,  Monterey,  California. 
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The  period  selected  (March  and  April  1972)  not  only  af- 
forded a  real-time  basis  for  investigation  but  it  is  also 
the  time  of  year  when  large,  well-defined  systems  charac- 
terize the  North  Pacific  Ocean.   These  systems  readily  lend 
themselves  to  the  modification  techniques  evaluated  during 
this  study . 
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IV.   PROCEDURES 

The  research  objectives  of  the  procedures  described 
below  were  to  extend  Bittner's  work  [1]  and  to  outline  an 
approach  which  would  be  practical  for  Fleet  use.   At  the 
time  when  the  FNWC  0000  GMT  numerical  analyses  are  received 
by  the  Weather  Centrals  and  Facilities,  the  prognostic 
fields  are  not  yet  available.   Therefore,  to  expedite  24- 
hour  cloud  forecasts  from  0000  GMT  it  would  be  best  to  use 
the  36-hour  prognostic  height  fields,  as  based  on  the  pre- 
vious 12-hour  analysis,  and  the  current  0000  GMT  analysis. 
Bittner  indicated  his  approach  is  good  for  24-  and  36-hour 
forecasts.   Since  he  effectively  assumes  that  there  is  no 
energy  transformation  occurring  between  the  SD  and  SR 
fields  for  24  to  36  hours,  then  a  change  in  an  SD  pattern 
may  be  propagated  forward  from  any  time  step  within  this 
forecast  period.   That  is,  a  change  in  the  SD  pattern  at 
plus  12  hours  from  the  time  of  the  forecast  may  be  carried 
forward  to  plus  36  hours.   Therefore,  in  these  procedures, 
changes  noted  in  the  0000  GMT  analysis  of  the  SD  patterns 
were  incorporated  in  the  prognostic  SD  patterns  which 
verified  24  hours  later,  but  which  were  made  from  the 
previous  1200  GMT  analysis. 

Another  point  which  may  be  mentioned  in  this  regard  is 
that  within  a  sparse  upper-air  data  area  such  as  the  North 
Pacific  Ocean,  the  difference  between  the  12-hour  prognosis 
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of  the  500-mb  level  and  the  verifying  500-mb  analysis  may 
be  quite  small. 

The  procedures  in  this  study  deal  almost  exclusively 
with  FNWC ' s  scale  separation  products  reproduced  on  a 
1:60,000,000  polar  st ereographic  projection.   This  required 
one  scale  change — that  of  manually  transferring  the  cloud 
features  of  the  ATS-1  photograph  to  this  scale.   FNWC  is 
currently  testing  a  program  that  will  convert  their  pro- 
ducts to  the  scale  and  projection  of  the  ATS-1  photograph. 
When  this  is  perfected,  all  future  studies  can  be  carried 
out  on  the  ATS-1  scale  to  eliminate  any  errors  in  the 
manual  transfer  process. 

In  general,  this  study  consisted  of  determining  a  modi- 
fied SD  field  and  projecting  these  modifications  on  the 
prognostic  SD  field.   Then  by  using  three  advecting  fields 
(FNWC's  Z  and  SR  PROGS,  and  a  hand-modified  FNWC  Z  PROG, 
called  Z'PROG),  three  cloud  cover  forecasts  were  determined 
for  verification.   The  detailed  procedures  are  as  follows: 

1)  Each  day  the  ATS-1  photograph  closest  in  time  to 
0000  GMT  was  selected  for  use.   The  major  cloud  features 
were  traced  onto  a  matching  grid  (same  scale  and  subpoint 
as  the  photograph),  carefully  delineating  as  covered  those 
areas  with  greater  than  50%  cloud  cover  and  as  open  those 
areas  with  less  than  50%  coverage. 

2)  The  cloud  patterns  were  then  transferred  by  hand  to 
a  polar-stereographic  projection  of  1:60,000,000  scale. 
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This  was  a  difficult  and  tedious  procedure  at  first  but 
after  some  practice,  the  time  was  decreased  to  less  than 
15  minutes  per  map. 

3)   FNWC ' s  SD  and  SR  analyses  were  put  in  register 
with  the  transferred  cloud  patterns.   In  areas  of  major 
cloud  features  such  as  fronts,  cyclones,  etc.,  the  SD 
patterns  and  isolines  were  modified  in  accordance  with 
the  rules  specified  by  Bittner  [1].   These  rules  are  listed 
in  Appendix  B  for  convenience  to  the  reader.   Further,  areas 
of  greater  than  50%  cloud  cover  were  assumed  to  be  associa- 
ted with  positive  vertical  motion  as  indicated  by  positive 
vorticity  advection  using  analyzed  SR  winds  as  the  advecting 
field.   The  SD  lines  were  redrawn  to  indicate  the  required 
advection  in  order  to  support  the  observed  cloud  patterns. 
Since  the  zero  lines  are  not  produced  on  the  SD  charts, 
these  lines  were  drawn  by  hand.   The  result  of  the  above  w.s 
a  modified  SD  chart,  hereafter  referred  to  as  the  SD' 
analysis.   The  SR  fields  were  not  modified. 

A)   Subjectively,  the  same  adjustments  that  were  made 
to  the  SD  analysis  were  then  applied  to  the  36-hour  SD 
prognosis,  hereafter  designated  as  the  SD '  prognosis.   Since 
the  location  of  zero-valued  SD  isolines  are  somewhat  arbi- 
trary, care  must  be  used  to  surround  the  same  centers  with 
the  same  zero  line  on  both  the  analysis  and  the  prognosis. 

5)   The  SD'  prognosis  was  then  graphically  added  to 
the  SR  prognosis  to  produce  a  modified  500-mb  height  prog- 
nosis, hereafter  indicated  as  a  Z'  prognosis. 
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6)  The  Z'  prognosis  was  put  in  register  with  the  SD' 
prognosis  and  used  as  the  advecting  field.   Clouds  were 
forecast  in  all  areas  of  negative  SD  advection  (positive 
relative  vorticity  advection).   Clear  skies  (less  than  50% 
covered)  were  forecast  for  areas  of  positive  SD  advection 
(negative  relative  vorticity  advection).   Note  that  clouds 
(greater  than  50%  covered)  were  assumed  to  be  located  only 
in  areas  of  positive  relative  vorticity  advection  in  this 
study.   This  is  a  slight  deviation  from  Bittner's  forecast- 
ing procedures  for  clouds  as  noted  below. 

7)  To  produce  the  second  cloud  cover  forecast,  steps 
1)  through  4)  are  the  same  but  step  5)is  omitted  since  the 
SR  prognosis  was  used  as  the  advection  field  for  the  pur- 
pose of  step  6) . 

8)  To  construct  the  third  cloud  cover  forecast,  steps 
l)through  4)  are  the  same  but  step  5)  is  omitted  since  the 
unmodified  Z  prognosis  was  used  as  the  advecting  field 
(Bittner's  approach  [1])  for  the  purpose  of  step  6). 

9)  Steps  1)  and  2)  above  were  then  applied  to  the 
next  day's  ATS-1  photograph  (current  plus  24  hours)  to 
produce  a  verifying  nephanaly sis . 

10)   The  verifying  nephanalysis  was  separately  put  in 
register  with  the  three  forecasts.  Each  of  the  351  grid 
points  for  each  of  the  three  forecasts  was  checked  and 
compared  with  the  corresponding  grid  point  of  the  verify- 
ing nephanalysis.   Grid  points  which  were  observed  clear 
but  forecast  cloudy  were  denoted  by  a  "2".   Those  which 
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were  both  forecast  and  observed  cloudy  were  indicated  with 
a  "0".   Next,  the  points  which  were  observed  cloudy  but 
forecast  clear  were  denoted  by  a  "1".   Finally  the  grid 
points  which  were  not  covered  by  clouds  on  the  verifying 
nephanalysis  nor  forecast  to  be  cloudy  were  represented  by 
a  "3".   At  this  point  the  technique  vrhich  gave  the  best 
results  could  be  readily  determined  by  totaling  the  number 
of  correct  forecasts  ("0"  and  "3")  for  each  forecast  day. 
Tables  I,  II  and  III  give  the  raw  scores  for  the  advecting 
fields  of  the  Z'  PROG,  SR  PROG  and  Z  PROG,  respectively. 
The  sixth  column  represents  the  number  of  correct  forecasts 
while  the  seventh  column  denotes  the  number  of  incorrect 
f orecas  t s . 

11)   Finally,  the  National  Meteorological  Center's 
threat  score  was  computed  for  each  of  the  forecasts.   An 
explanation  of  the  threat  score  is  presented  in  Section  VI. 

As  indicated  earlier,  the.  manner  in  which  the  clouds 
were  determined  by  this  author  is  not  strictly  the  same  a? 
that  used  by  Bittner  [1].   In  addition,  he  applied  the  rules 
of  Appendix  B  in  reverse  to  indicate  some  geometric  detail 
in  the  cloud  patterns,  as  spiral  clouds,  vortices,  etc., 
whenever  they  were  deemed  necessary  as  determined  by  the 
SD'  prognosis.   He  also  used  various  crossing-angle  rules 
as  a  guide  for  forecast  cloud  tops,  areas  and  amounts  of 
precipitation,  etc.   Since  this  study  was  conducted  to  try 
and  improve  Bittner's  basic  premise,  these  refinements  were 
not  included . 
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V .   CASE  STUDY  OF  MODIFICATIONS  TO  ANALYSES 

The  synoptic  time  of  0000  GMT  10  April  1972  and  its 
associated  analyses  and  prognoses  were  randomly  selected 
as  an  example  of  the  procedures  presented  in  the  preced- 
ing section.   Figures  1  through  18  graphically  illustrate 
the  following  discussion. 

Figure  1  is  the  ATS-1  photograph  of  2343  GMT  9  April 
1972  which  was  selected  to  be  representative  of  the  cloud 
cover  for  0000  GMT  10  April  1972.   The  major  cloud  fea- 
tures (areas  of  >_  50%  cloud  cover)  are  outlined  for 
clarity. 

In  general,  SD  patterns  moving  into  the  western  sector 
of  the  verification  grid  were  modified  very  little  because 
the  cloud  patterns  are  not  well-defined  near  170°E  in  the 
ATS-1  pictures.   This  sometimes  led  to  a  large  error  in 
the  verification  of  cloud  forecasts  west  of  180°.   There 
was  usually  good  definition  of  the  cloud  features  and, 
therefore,  accurate  modifications  of  the  SD  patterns  in 
the  extreme  eastern  section  of  the  grid.   These  eastern 
SD  features  and  associated  cloud  features  sometimes  moved 
out  of  the  area  of  the  verification  grid  during  the  time- 
span  of  the  prognosis.   Some  of  the  centers  of  the  SD 
analyses  in  this  near  North  American  coastal  area  were 
extensively  modified.   For  those  systems  which  did  not 
appear  within  the  verification  area  at  the  verifying  time 
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of  the  prognosis,  the  overall  goodness  of  the  associated 
modifications  and  the  resultant  cloud  forecasts  are  un- 
known . 

Figure  2  shows  a  nephanalysis  (areas  of  >_  50%  cloud 
cover  as  derived  from  Figure  1,  on  a  grid  of  the  same 
scale . 

Figure  3  is  the  nephanalysis  after  it  has  been  trans- 
ferred to  a  polar  s tereographic  projection  with  the  SD 
and  SR  fields  superimposed.   The  SD  pattern  is  illustrated 
with  dashed  lines  and  the  SR  field  with  solid  lines.   Note 
the  strong  SD  low  center  at  43°N  and  135°W.   Such  a  low 
value  for  the  SD  center  suggests  a  cyclone  at  500  mb . 
Since  this  center  is  embedded  in  an  SR  trough,  as  indicated 
by  the  solid  lines,  there  is  high  probability  of  a  closed 
center  at  500  mb .   However,  the  cloud  patterns  of  Figure  1 
indicate  only  slight  cyclonic  circulation  with  convective 
clouds  in  the  area.   This  suggests  decreasing  the  positive 
vorticity  advection  (negative  SD  advection)  considerably. 

Figure  4  shows  the  new  or  modified  SD  field.    Note 
the  drastic  changes  made  to  the  SD  zero  line.   In  the  lower 
left  portion  of  the  figure,  the  zero  line  was  relocated  to 
the  poleward  edge  of  the  frontal  cloud  band.   (Rule  l.a 
Appendix  B) .   In  general,  SD  troughs  and  ridges  were  re- 
drawn to  coincide  with  major  cloud  boundaries.  Note  the 
considerable  filling  of  the  SD  low  center  at  43°N,  135°W 
to  decrease  the  negative  SD  advection.  The  center  at  33°N 
180°  was  moved  slightly  northward  to  maintain  the  correct 
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perspective  with  the  circulation  center  on  the  photograph. 
(See  rule  2.e  Appendix  B.) 

Figure  5  is  the  modified  SD  analysis  of  Figure  A 
(dashed  lines)  superimposed  on  the  original  SD  field  of 
Figure  3  (solid  lines).   The  amount  and  direction  of  move- 
ment of  all  centers  and  isolines  were  noted  and  used  in 
the  modification  of  the  SD  prognosis. 

Figure  6  represents  the  36-hour  SD  prognosis  and  its 
modification  (SD1  PROG).   As  such,  it  represents  the 
changes  made  to  the  SD  analysis  24  hours  earlier,  as  shown 
in  Figure  5.   As  can  be  seen  from  the  figure,  this  modifi- 
cation is  somewhat  subjective. 

Figure  7  shows  the  graphical  addition  of  the  modified 
SD  prognosis  (SD'  PROG)  to  the  SR  prognosis  in  order  to 
produce  a  new  height  prognosis  (Z'  PROG).   This  is  the 
only  additional  step  required  over  Bittner's  method. 

Figure  8  is  the  Z'  prognosis  of  Figure  7  superimposed 
on  the  SD'  prognosis  of  Figure  6.   Since  the  negative  SD 
values  correspond  to  positive  values  of  vorticity,  as 
shown  in  Appendix  A,  the  areas  of  positive  vorticity  advec- 
tion  (PVA)  are  readily  apparent .   These  areas  are  carefully 
delineated  and  are  represented  in  the  figure  as  clouds. 
Note  that  major  cloud  boundaries  are  determined  by  the  SD 
troughs  and  ridges  as  defined  earlier  relative  to  vorticity 
advection  and  as  stressed  by  Bittner.   As  mentioned  in  the 
previous  section,  Bittner  inserted  cloud-pattern  detail 
(i.e.  cloud  vortices,  developing  cyclones,  etc.)  at  this 


28 


point.   These  modifications  could  be  done  by  a  knowledge- 
able synoptic  meteorologist  resulting  in  better  continuity 
of  cloud  patterns  and  increased  accuracy  of  cloud  fore- 
casts . 

Figures  9  and  10  are  the  same  as  Figure  8  except  that 
the  vor t icity-advect ion  areas  are  determined  by  using  the 
SR  unmodified  Z  prognoses  as  the  advection  fields,  respec- 
tively.  See  procedures  7  and  8  of  the  previous  section. 

Figure  11  is  the  verifying  ATS-1  photograph  with  areas 
of  >_  50%  cloud  cover  outlined.   This  picture  also  shows 
very  little  cyclonic  motion  in  the  northeast  corner,  justi- 
fying the  decreased  PVA  of  Figure  4. 

Figure  12  shows  the  boundaries  of  the  >_  50%  cloud  cover 
from  the  photograph  of  Figure  11. 

Figure  13  illustrates  the  verifying  clouds  as  trans- 
ferred from  Figure  12  and  the  PVA  areas  derived  in  Figure 
8  from  the  Z1  prognosis.   As  mentioned  earlier,  the  date 
for  this  case  study  was  chosen  at  random;  it  just  happens 
to  be  a  below-average  date.   See  Tables  I,  II  and  III  for 
comparisons  of  raw  values  for  0000  GMT  11  April  1972. 

Figure  14  represents  the  accuracy  of  the  forecast. 
The  351  grid  points  of  Figure  13  (every  2-1/2°  latitude- 
longitude  spacing)  were  examined  and  assigned  to  one  of 
four  categories.   The  grid  points  that  were  forecast  to 
be  cloudy  but  were  observed  to  be  clear  ("2")  totaled  65. 
Those  that  were  both  forecast  and  observed  cloudy  ("0") 
amounted  to  141.   Observed  cloudy  but  forecast  clear  ("1") 
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were  101  points.   And,  finally,  those  grid  points  which 
were  both  forecast  and  observed  to  be  clear  ("3")  totaled 
A3.   This  amounted  to  184  correct  ("0  +  3")  and  166  ("1  + 

2")  incorrect  forecasts,  which  yields  a  combined  threat 

2 
score   of  .67.   See  discussion  in  Section  VI  and  scores 

in  Table  IV. 

Figures  15  and  16  show  the  verifying  clouds  superim- 
posed on  the  PVA  area  derived  from  the  SR  prognosis  and 
the  scoring  grid,  respectively.   The  number  of  grid  points 
with  correct  forecasts  ("0  +  3")  amounted  to  170  while  the 
incorrect  total  ("1  +  2")  was  181.   The  threat  score  for 
this  method  was  .59. 

In  general,  the  scores  for  all  three  methods  in  this 
case  study  were  below  their  respective  averages.   See 
Table  IV  for  the  averages.   However,  the  relative  differ- 
ences between  the  three  methods  are  of  the  same  order  as 
the  overall  averages  of  each.   That  is  the  SR  method  of 
determining  the  PVA  areas  yielded  approximately  the  same 
score  as  Bittner's  technique  (unmodified  Z  field)  for 
determining  PVA  areas  while  the  Z'  (modified  Z  field) 
m  ethod  gave  better  than  a  10%  increase  in  the  value  of  the 
combined  threat  score  which  was  used  to  evaluate  the  cloud 
cover  forecasts. 


Threat  score  will  be  explained  in  detail  in  the  following 
sect  ion . 
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VI.   VERIFICATION  OF  RESULTS 


The  National  Meteorological  Center's  threat  score  sys- 
tem was  used  to  measure  the  goodness  of  the  forecasts. 
The  score  depends  on  the  correct  forecasting  of  cloudy 

(T_  )  and  clear  areas  (T.) .   They  are  calculated  as  follows 
w  U 


0 


• 


0+1  +  2 


and    T„  ■ 


0    3+1+2 


where : 


0  =  Those  points  both  forecast  and  observed  as 

cloudy . 

1  =  Those  points  forecast  clear  but  observed  as 

cloudy  . 

2  •=  Those  points  forecast  cloudy  but  observed  as 

clear . 

3  =  Those  points  both  forecast  and  observed  as 

clear . 


Using  this  approach,  a  perfect  forecast  of  clouds  re- 
sults in  a  score  of  T    =  1.0  and  a  perfect  forecast  of 
clear  skies  results  in  a  score  of  T   =  1.0.   Combining 
these  two  perfect  threat  scores  gives  a  total  of  2.0  for 
an  absolutely  correct  forecast  over  an  area  as  follows: 

T^  +  T.  =  1.0  +  1.0  =  2.0 
V-  0 

The  results  of  the  three  forecasts  prepared  for  each 

day  are  listed  in  Table  IV.   These  scores  were  obtained 

using  the  above  technique.   A  cloudy  forecast  score  and 

a  clear  forecast  score  followed  by  a  combined  score  are 
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presented  for  each  day  of  the  study  period.   The  overall 
averages  of  each  method  are  listed  at  the  bottom  of  the 
table. 

The  results  of  Bittner's  study  of  approximately  the 
same  area  of  the  North  Pacific  Ocean,  but  for  different 
dates,  are  presented  in  Table  V  for  comparison.   The  un- 
modified column  represents  threat  scores  obtained  without 
the  use  of  satellite  pictures  to  modify  either  the  initial 
or  prognostic  SD  fields.   The  modified  column  represents 
the  forecasts  obtained  after  the  initial  SD  fields  have 
been  modified  based  on  the  satellite  observed  cloud  pat- 
terns and  the  prognostic  SD  fields  have  been  modified 
based  on  the  modifications  applied  to  the  initial  SD 
patterns.   This  latter  column  represents  a  15%  increase 
in  cloud-cover  forecasting  accuracy  over  the  unmodified 
column.   Except  for  the  additional  refinement  steps  men- 
tioned earlier,  this  modified  column  is  equivalent  to  the 
Z  PROG  column  of  Table  IV.   The  additional  refinements  in 
forecasting  cloud  cover  used  by  Bittner  may  account  in 
large  part  for  the  relative  accuracy  differences  between 
the  two  corresponding  columns  of  Tables  IV  and  V.   How- 
ever, as  noted  earlier,  these  refinements  may  be  easily 
added  to  the  basic  approach  in  determining  the  PVA-associa- 
ted  cloud  cover. 

The  basic  determination  of  the  PVA  areas  by  the  three 
different  methods  used  in  this  study  are  clearly  comparable 
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The  results  in  Table  IV  indicate  that  the  Z  PROG  and  the 
SR  PROG  advection  fields  are  about  equal  in  their  results. 
The  Z'  PROG  advection  fields  indicate  about  a  10%  improve- 
ment over  the  approach  suggested  by  Bittner  [1].   Thus  it 
may  be  argued  that  if  Bittner's  approach  is  modified  to 
the  extent  of  including  the  Z'  (vice  the  Z)  prognosis  for 
the  advecting  field,  before  adding  the  required  and  obvious 
continuity  of  large-scale  cloud  patterns,  the  best  24-hour 
cloud-cover  fields  may  be  produced. 
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TABLE  IV.   A  comparison  of  threat  scores  of  Z'PROG,  SR 
PROG  and  Z  PROG  cloud-cover  forecasting  techniques. 


Verificat 

:ion 

DATE  /Time 

Z  'PROG 

SR  PROG 

Z  PROG 

1972 

00  GMT 

T 

• 

+  To 

=  T 

T 

• 

+  T 
0 

=  T 

T0 

+  To 

=  T 

3/10 

.95 

.51 

1.46 

.54 

.48 

1.02 

.49 

.46 

.95 

3/15 

.41 

.25 

.66 

.37 

.23 

.60 

.36 

.22 

.58 

3/16 

.51 

.43 

.94 

.52 

.41 

.93 

.50 

.42 

.92 

3/17 

.A3 

.40 

.83 

.44 

.40 

.84 

.45 

.43 

.88 

3/22 

.33 

.36 

.69 

.34 

.36 

.  70 

.30 

.37 

.67 

3/23 

.48 

.48 

.96 

.47 

.49 

.96 

.46 

.48 

.94 

3/24 

.54 

.28 

.82 

.47 

.25 

.72 

.50 

.27 

.77 

3/25 

.53 

.34 

.87 

.50 

.33 

.83 

.49 

.32 

.81 

A/1 

.67 

.45 

1.12 

.55 

.43 

.98 

.53 

.36 

.89 

A/A 

.41 

.22 

.63 

.51 

.22 

.73 

.52 

.18 

.70 

A/5 

.61 

.14 

.75 

.53 

.12 

.65 

.60 

.15 

.75 

A/6 

.57 

.14 

.71 

.49 

.14 

.63 

.55 

.16 

.71 

A/7 

.76 

.22 

.98 

.72 

.18 

.90 

.70 

.18 

.88 

A/8 

.67 

.29 

.96 

.63 

.23 

.86 

.62 

.34 

.96 

A/9 

.63 

.44 

1.07 

.63 

.43 

1.06 

.59 

.38 

.97 

A/10 

.59 

.29 

.88 

.  75 

.21 

.96 

.47 

.18 

.65 

A/11 

.46 

.21 

.67 

.43 

.16 

.59 

.41 

.18 

.59 

A/12 

.89 

.15 

1.04 

.53 

.16 

.69 

.53 

.19 

.72 

A/13 

.44 

.32 

.76 

.43 

.34 

.77 

.38 

.32 

.  70 

A/14 

.27 

.24 

.51 

.26 

.24 

.50 

.34 

.25 

.59 

Average 

.56 

.31 

.87 

.51 

.29 

.80 

.49 

.29 

.  78 
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TABLE  V.   The  results  of  Bittner's  study  [1] 


DATE 


UNMODIFIED 


MODIFIED 


1969 


+  To   = 


+  T„  =  T 


10/3 

.34 

.66 

1.00 

10/9 

.30 

.63 

.93 

10/10 

.29 

.49 

.78 

10/13 

.24 

.58 

.82 

10/14 

.14 

.60 

.74 

10/15 

.28 

.59 

.87 

10/16 

.24 

.57 

.81 

10/17 

.35 

.68 

1.03 

10/18 

.24 

.57 

.81 

10/19 

.45 

.79 

1.24 

10/20 

.27 

.62 

.89 

10/21 

.32 

.65 

.97 

10/24 

.26 

.43 

.69 

10/25 

.38 

.61 

.99 

10/27 

.25 

.60 

.85 

10/30 

.26 

.35 

.61 

11/4 

.32 

.57 

.89 

.22 

.61 

.83 

.  33 

.  73 

1.06 

.34 

.66 

1.00 

.32 

.66 

.98 

.23 

.63 

.86 

.32 

.68 

1  .00 

.34 

.74 

1.08 

.26 

.70 

.96 

.36 

.80 

1.16 

.40 

.84 

1.24 

.31 

.68 

.99 

.34 

.71 

1.05 

.25 

.68 

.93 

.27 

.60 

.87 

.43 

.78 

1.21 

.40 

.63 

1.03 

.28 

.57 

.85 

Average    .29   .59 


.88 


.32   .69   1.01 
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VII .   CONCLUSIONS 

This  study  indicates  there  is  merit  in  using  the  satel- 
lite-observed cloud  patterns  to  modify  the  position  and 
orientation  of  500-mb  relative  vorticity  patterns  (the  SD 
fields  of  FNWC)  in  order  to  improve  24-hour  forecasts  of 
cloud  cover  over  oceanic  areas.   Bittner  [1],  in  an  earlier 
study,  came  to  the  same  conclusion. 

It  appears  that  Bittner's  basic  approach  can  be  improved 
by  using  both  modified  prognostic  relative  vorticity  pat- 
terns and  modified  prognositic  height  fields  at  500  mb . 
However,  the  refinements  of  cloud-pattern  continuity  de- 
scribed by  Bittner  should  be  included  in  the  final  24-hour 
cloud  pattern  forecast  for  best  results. 
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Figure  1.   2343  GMT  09  April  1972  ATS-1  satellite  photo- 
graph (Subpoint:  2°S).   Area  of  >^  50%  cloud  cover  bounded 
by  3    . 


Figure  2.   Nephanalysis  of  ATS-1  satellite  photograph 
(Subpoint:  2°S),  2343  GMT  09  April  1972.   Area  of  >  50% 
cloud  cover  bounded  by  £?  . 
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Figure  3.   FNWC  500-mb  SD  ( )  and  SR  ( )  analyses,  0000 

GMT  10  April  1972,  superimposed  on  the  ATS-1  nephanalysis 
of  >_  50%  cloud  cover  (2343  GMT  09  April  1972).  SD  interval: 
30  m;  zero  SD  line  labeled. 
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Figure  4.   FNWC  500-mb  SR  ( )  and  modified  SD  ( ) 

analyses,  0000  GMT  10  April  1972,  superimposed  on  the 
ATS-1  nephanalysis  of  >_  50%  cloud  cover  (2343  GMT  09 
April  1972).  SD  interval:  30  m;  zero  SD  line  labeled. 
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Figure  5.   500-mb  FNWC  SD  ( )  and  modified  FNWC  SD  ( ) 

analyses,  0000  GMT  10  April  1972:  SD  interval:  30  m;  zero 
SD  line  labeled. 


Figure  6.   36-hour  500-mb  FNWC  SD  ( )  and  modified  FNWC 

SD  ( )  prognosis,  verifying  at  0000  GMT  11  April  1972: 

SD  interval:  30  m;  zero  SD  line  labeled. 
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Figure  7.   Modified  FNWC  Z  prognosis  ( )  as  derived  from 

the  FNWC  SR  (-•-•)  and  modified  SI)  ( )  prognosis,  verify- 
ing at  0000  GMT  11  April  1972:  SD  interval:  30  m;  zero  SD 
line  labeled. 


Figure  8.   Modified  FNWC  500-mb  Z  ( )  and  SD  ( )  prog- 
noses with  resulting  areas  of  positive  vorticity  advection 
(  Q     ),  verifying  at  0000  GMT  11  April  1972:  SD  interval: 
30  m;  zero  SD  line  labeled. 
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Figure  9.   FNWC  500-mb  SR  ( )  and  modified  SD  ( )  prog- 
noses with  resulting  areas  of  positive  vorticity  advection 
(  O     ) ,  verifying  at  0000  GMT  11  April  1972:  SD  interval: 
30  m;  zero  SD  line  labeled. 


Figure  10.   FNWC  500-mb  Z  ( )  and  modified  SD  ( )  prog' 

noses  with  resulting  areas  of  positive  vorticity  advection 
(  0     ).  verifying  at  0000  GMT  11  April  1972:  SD  interval: 
30  m;  zero  SD  line  labeled. 
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Figure  11.   2358  GMT  10  April  1972  ATS-1  satellite  photo- 
graph (Subpoint:  2°S).   Area  of  >_  50%  cover  bounded  by  ££. 


Figure  12.   Nephanalysis  of  ATS-1  satellite  photograph 
(Subpoint:  2°S),  2358  GMT  10  April  1972.   Area  of  >  50% 
cover  bounded  by  £3 
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Figure  13.   Nephanalysis  of  >_  50%  cloud  cover  from  ATS-1 
photograph  (2358  10  April  1972)  (  £3  )  ,  and  area  of  posi- 
tive vorticity  advection,  (:a£V:  )  derived  from  the  modified 
36-hour  FNWC  500-mb  Z  prognosis  verifying  at  0000  GMT 
11  April  1972. 
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Figure  14.   Scoring  grid  for  the  cloud  forecast  based  on  the 
modified  36-hour  FNWC  500-mb  Z  prognosis,  verifying  at  0000 
GMT  11  April  1972.   Legend:  0-Clouds:  forecast  and  observed- 
1A1  points;  1-Clouds:  observed  but  not  forecast-101  points; 
2-Clouds:  forecast  but  not  observed-66  points;  3-Clear:  fore- 
cast and  observed-43  points. 
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Figure  15.   Nephanalysis  of  >^  50%  cloud  cover  from  ATS-1 
photograph  (2358  GMT  10  April  1972)  ( Q       ),  and  area  of 
positive  vorticity  advection  derived  from  36-hour  FNWC 
500-mb  SR  prognosis  verifying  at  0000  GMT  11  April  1972. 
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Figure  16.   Scoring  grid  for  the  cloud  forecast  based  on  the 
modified  36-hour  FNWC  500-mb  SR  prognosis,  verifying  at  0000 
GMT  11  April  1972.   Legend:  0-Clouds:  forecast  and  observed- 
135  points;  l«»Clouds:  observed  but  not  forecast-105  points; 
2-Clouds:  forecast  but  not  observed-76  points;  3-Clear:  fore- 
cast and  observed-35  points. 
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Figure  17.   Nephanalysis  of  >^  50%  cloud  cover  from  ATS-1 
photograph  (2358  GMT  10  April  1972)  ( <?   ),  and  area  of 
positive  vorticity  advection  derived  from  the  36-hour 
FNWC  500-inb  Z  prognosis  verifying  at  0000  GMT  11  April  1972. 
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Figure  18.   Scoring  grid  for  the  cloud  forecast  based  on  the 
modified  36-hour  FNWC  500-mb  Z  prognosis,  verifying  at  0000 
GMT  11  April  1972.   Legend:  O-Clouds:  forecast  and  observed- 
127  points;  1-Clouds:  observed  but  not  forecast-110  points; 
2-Clouds:  forecast  but  not  observed-74  points;  3— Clear:  fore 
cast  and  observed-40  points. 


45 


B 


foOnhSO     -£,'(/ 


Oo  AS       ltht 
'  Cen-t<   fS 


Figure  19.   The  relationship  between  cloud  patterns  and  the 
500-mb  SD  field,  as  presented  by  Bittner  [1]  and  illustrated 
by  Zeigler  [8]:  A.  Surface  wave  vs.  SD  field.  B.  Occluded 
system  with  500-mb  flow  and  SD  field.  C.  Mature  cyclone  vs. 
SD  field.  D.  Dissipating  cyclone  vs.  SD  field.  E.  Polar  jet 
axis  and  open-celled  cumulus  vs.  zero  SD  contour.  F.  SD  flow 
vs.  cloud-associated  PVA  MAX. 
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APPENDIX  A 

FNWC'S  500-MB  SD  FIELD  AS  AN  APPROXIMATION  TO 
THE  RELATIVE  VORTICITY  FIELD 


The  definitions  and  discussion  below  follow   closely 
that  in  [1]  . 

The  500-mb  disturbance  component  is  defined  as 

SD500  E  Z500-  SR500  =  -"2  '"         V2z500da         (1> 

O 

where : 

Z     =  500-mb  height  field 

SR  nn  =  500-mb  residual  (mostly  wave  numbers  1,2,3,4, 
500    5  and  6) 

SD-flfl  ■  500-mb  shortwave  field  (mostly  wave  numbers 
7  and  above) 

a  =  degree-of -smoothing  parameter 

SS  =  numerical  grid  length 

The  right  hand  side  of  this  equation  is  a  Laplacian  of  the 

height  field  over  the  smoothing  interval  a  . 

Geostrophic  relative  vorticity  is  defined  as: 

5  =  gf"1  V2Z50Q  (2) 

where:  g  =  acceleration  due  to  gravity 

f  =  Coriolis  parameter 
By  substituting  Equation  (2)  into  Equation  (1),  the  follow- 
ing relationship  is  obtained: 

o    _1    aSD 
SD5Q0  =  -6SZ  fg  X   /     Ua  (3) 

o 
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Equation  (3)  shows  that  the  SD  field  is  about  the  same  as 
the  geostrophic  relative  vorticity  field  (with  a  sign  re- 
versal) over  the  smoothing  interval  used  in  the  analysis. 
Thus,  for  all  practical  purposes,  they  can  be  considered 
analogous  for  both  analyses  and  prognoses. 
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APPENDIX  B 

SUMMARY  OF  RULES  RELATING  CLOUD 
PATTERNS  TO  500-MB  SD  ANALYSES 


Following  is  a  summary  of  rules  relating  cloud  patterns 
to  the  500-mb  analyses  [1]. 

1)  The  SD  zero  line  should  be  drawn:   (See  Figure  19  for 
a  graphical  illustration  of  some  of  these  rules.) 

a)  Along  the  poleward  edge  of  the  cloud  band  associated 
with  the  polar  front. 

b)  Along  or  just  equatorward  of  the  cloud  shadow  mark- 
ing the  polar  jet  axis;  along  or  just  inside  the 
poleward  limit  of  transverse  lines  marking  the  sub- 
tropical-jet axis. 

c)  Along  the  line  between  two  major  vorticies  which 
are  connected  by  a  continuous  frontal  band. 

2)  An  SD  low  center  should  be  located: 

a)  Five  to  seven  degrees  upstream  from  an  open  frontal 
wave,  or  at  the  rear  edge  of  the  head  of  a  cloud- 
identified  positive  vorticity  advection  maximum 
(PVA  MAX) . 

b)  Five  degrees  upstream  from  the  center  of  a  hook- 
shaped  cloud  pattern  (rapidly  developing  cyclone). 

c)  Within  three  degrees  or  less  from  the  center  of  a 
complete  spiral  cloud  pattern  (mature  cyclone). 

d)  Directly  over  or  slightly  to  the  east  of  the  center 
of  an  open  vortex  (dissipating  cyclone). 
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e)   Over  the  center  of  a  strong  vorticity  center  in 
the  cold  air . 
3)   An  SD  trough  should  be  drawn: 

a)  Through  the  point  where  a  frontal  cloud  band 
forks  or  fractures. 

b)  Five  to  seven  degrees  upstream  from  a  weak  frontal 
wave  (where  a  closed  SD  low  is  not  justified). 

A)   An  SD  ridge  should  be  drawn: 

a)  In  an  occlusion,  to  coincide  with  the  downstream 
major  cloud  boundary  or  the  change  from  C  (cumuli- 
form/stratiform)  to  MCO  (s trat if orm/cirrif orm) . 

b)  With  anticyclonic  curvature  near  the  equatorward 
edge  of  a  frontal  cloud  band. 

c)  With  anticyclonic  curvature  over  the  peak  of  a 
frontal  wave  pattern. 

5)  The  SD  gradient  should  be  maximized: 

a)  Along  the  rear  edge  of  a  frontal  cloud  band. 

b)  In  the  dry  tongue  of  a  hook-shaped  vortex. 

c)  At  and  poleward  of  the  peak  of  a  frontal  wave. 

6)  SD  lines  should  be  drawn: 

a)  Parallel  to  lines  of  cumulus  (CELLS  OPEN)  in  the 
cold-air  sector  of  a  cyclonic  disturbance. 

b)  Parallel  to  cloud  lines  and  bands  in  a  dissipating 
vortex . 
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